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RESO SUCHOn: 

As a part of the problem of developing numerical wave 
forecasting procedures for the North Atlantic Ocean, selected se- 
quences of the weather maps for the North Atlantic for which wave 
data were known to be available were studied in detail for the five 
year period beginning in April 1955 and ending in March 1960. 
Certain dates and times of observations were selected for a variety 
of reasons for study. For these dates and times, the National Insti- 
tute of Oceanography provided copies of the wave records that were 
obtained by the OWS Weather Explorer and by the OWS Weather Reporter. 

In total, about 800 wave records were provided, anda 
complete spectral analysis is planned for about 400 of these records. 

This report is the first of a series of reports to present 


in tabular and graphical form the results of these analyses. The total 


number of spectra given is 114. 


The original wave records varied in length, but almost all 
of this first set were 15 minutes long. The crest to trough heights 
of the highest waves in a particular record (uncorrected for cali- 
bration effects) varied from a few feet to more than 60 feet in the 
complete set of records. Bounds were set on each record just above 
the highest wave crest and just below the lowest wave trough, and 
the records were read to an accuracy of one part in a thousand 
(nominally) over this range at an interval of 1.5 seconds throughout 


the record. Thus a 15 minute record was reduced to a time series 
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of 600 points. Where gaps or irregularities occurred, the records 
were smoothed by hand as accurately as possible. 

The time series of 600 points was then analyzed on the 
CDC 1604 so as to estimate the energy spectrum of the waves at 60 
points over the frequency range from zero to 0.333 cycles per second 
by means of the procedures given by Tukey (1949) as explained in 
detail by Blackman and Tukey (1958). The smoothing operation that 
was used to go from L to U inthe equations of Blackman and Tukey 


was 


(1) UL = 0.251, 1 + 0.50L, + 0.25L,, 


1 
with suitable corrections at the ends of the range. 

The spectral estimates so obtained still had to be cor- 
rected for the response of the shipborne wave recorder (Tucker, 
1956) and for the introduction of noise in both the original record 
and in the digitization procedure. The calibration of the shipborne 
recorder depends on the ship, and the calibration curves were pro- 
vided by Mr. D.E. Cartwright for this purpose. The calibration curve 
for the Weather Explorer is given by Table 1. The calibration curve 
for the Weather Reporter is given by Table 2. 

As in another investigation (Bretschneider, Crutcher, et 
al (in press)), it was fiend eR the application of the above cali- 
bration curves to the spectra that were originally computed resulted 
first in a decrease and then a rapid increase in the spectra at high 
frequencies due to the presence of noise and other irregularities 
(possibly from nonlinear effects in the original wave records) at 


the high frequency end of the spectrum. To eliminate this effect, 
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the last part of the spectrum was smoothed by a three point running 
weighted mean (0.25, 0.50, and 0.25) and then the last ten values 
were averaged.” This average was treated as white noise and sub- 
tracted from all spectral estimates. When the reduced values were 
multiplied by the appropriate calibration curves, the usual result 
was a fairly smooth spectrum that decreased regularly toward zero 
values at high frequency. By such a procedure some of the spectral 
values at high frequency will be negative. These values were auto- 
matically set equal to zero in the rest of the computations. 

Inadvertently, Table 1 was applied where Table 2 ought 
to have been applied to some of the spectra. The result was ex- 
ceptionally high values at high frequencies. Thus the calibrations 
given in these two tables do distinguish between the characteristics 
of the two ships as, when done correctly, reasonable results were 
obtained. 

Even with these corrections, there were a few spectra that 
still became exceptionally large for frequencies greater than about 
0.25 cycles per second. This behavior was apparently caused by the 
original quality of the record and not by the digitization procedure. 
These spectra were further modified by arbitrarily setting the cali- 
bration curve equal to one above a certain frequency that was 
selected by inspection of each spectrum. 

The result of such a sequence of computations should yield 
fairly reliable spectral estimates for frequencies ranging from zero 
to 0.25 cycles per second, but the values at high frequencies should 


not be used to decide on any features of the high frequency end. 


*See p. 6 for details. 
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The spectral estimates that resulted from this sequence of 
operations were then processed further to obtain some additional 
useful information. Let Ue for h=0, 1, 2, --+,60, represent the 
spectral estimates (after subtraction of the noise and multiplication by 
the calibration for the shipborne recorder) in terms of the resolution 


of the variance of the wave record into frequency intervals. The 


following quantities were then also computed and tabulated with each 


spectrum. 

(2) CORR VAR = corrected variance = 72100 
== *,1/2 

(3) SIG HGT = Hi =E2OS (22U, ) 

4 AVERT = ¥ = (50, */=£20,*]" 

(4) a = PU aaa 


(5) TOTAL DF= Total degrees of freedom=10 [2 cavilal [2 OT 


(for 600 points, 60 lags; i.e., 20 degrees of freedom 
per spectral estimate) 
The confidence intervals on the corrected variance and on the 


significant height are given by 


Upper 95% on CORR VAR =(10"!/NTPF ) CoRR VAR 
(6) 
evens, on CORRI VAR = (lh. > yCoRR VAR 


and by 
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3 3 
in terms of the total degrees of freedom (TDF) to a high degree of 
accuracy since the total degrees of freedom are large. 

The corrected variance, the significant height, and the 
total degrees of freedom are relatively insensitive to changes in the 
noise level and in the high frequency behavior of the spectrum. 
However, the average period can properly be viewed with caution. 

The winds near the ship at the time of observation are 
also given to the nearest five knots as read directly from weather 
maps. These values are subject to later correction in terms of the 
logs of the weather ships. 

He lana voureptablesiand gs nepns) 

The body of this report consists of supplementary tables, 
of tables that give the appropriate results for each of the original 
wave records, and of graphs of each of the estimated spectra along 
with the confidence intervals on the spectra. 

The supplementary tables consist of Tables 1 through 4. 
Tables 1 and 2 have been described above. 

Table 3 gives either the on station position of the ship, 
A, I, J, or K, or the latitude and longitude of the ship if it is going 


on or off station. The speed and direction of the ship is given. 
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Position A corresponds to 62°N, 33°W. 
Position I corresponds to 59°N, 19°W. 
Position J corresponds to 52.5°N, 20°W. 
Position K corresponds to 45°N, 16°W. 


If the record was not 15 minutes long, less than 600 points were 


read. For these records, Table 4 gives the actual number of points 


used and the corrected total degrees of freedom. A correction to the 


upper and lower confidence limits, which would be quite small, would 


also be needed to be exact. 
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A tabulated spectrum can be interpreted as follows: 

la) Supplementary data for each spectrum consist of the date, 
hour, wind speed, total degrees of freedom, average 
period, significant height, corrected variance, noise 
level, and record number. Some tables give the confi- 
dence limits for the height according to equation (7). 

lb) In the first column, the spectral lag numbers (H) are given. 

2) In the second column (FRE) the frequency according to the 
equation f = H/180 (sec) is given. 

3) In the third column (UNIT = FT*), the spectrum as com- 
puted from the original data is given in units of (£t)”. 

4) Inthe fourth column (FILTERED), a smoothing operator 


for H> 40 is applied. It is actually 


Puy = 0.25U,,4 + 0.50U,, + 0.25F i) 
(where F = Filtered, and U = Unit) 
5) In the fifth column (LESS NOISE), the noise level shown at 


the top is subtracted from each estimate. 
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6) In the sixth column (CORR FT 2), the LESS NOISE 
column is multiplied by the calibration curve for the 
shipborne record according to either Table 1 or 
Table 2. If this column agrees with the previous column, 
at high frequencies, the calibration curve has been arbi- 
trarily set equal to one to avoid extreme values at high 
frequency. 

7) Inthe last two columns, the upper and lower 95% and 5% 


confidence bounds are shown. 
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The graph that accompanies the spectral tabulation shows the 
spectrum and the 95% and 5% confidence bounds. The scale is chosen 
so that the highest 95% confidence value is at the top of the graph and 
the vertical axis of the coordinate system shows the spectral values 
for that spectrum in units of (ects: The scales change with each 
spectrum and comparisons between spectra by means of the graphs 


should be made cautiously. 
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Table 1. Calibration factors for the Weather Explorer. 


1.0000 

1.0C00 1.0000 1.0000 1.0000 1.0006 1.6157 
1.3740 1.2452 1.1746 1.1399 1.1291 1.1343 
1.1547 1.1870 1.2304 1.2845 1.3504 1.4277 
1.5193 1.6241 1.7444 1.8828 2.0415 222243 
224349 2.26765 22-9523 322725 326414 4.0714 
4.5654 5.1490 528190 6.6136 1.25383 8.6336 


9.9169 11.4459 13.2691 15.4245 18.0095 21.1086 
24-8366 29.3522 34.8079 41.4485 49.5464 5924548 
71.5502 86.5947 105.1503 128.1186 156.7723 192.5202 

237-3987 293.8682 365.1736 455.5306 570.2699 716.8705 
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Table 2. Calibration factors for the Weather Reporter. 


1.0000 

1.0000 1.0000 1.0000 1.C000 1.0000 1.5755 
1.3277 1.1908 1.1099 1.0630 1.0375 1.0257 
1.0260 1.0350 1.0514 1.0633 1.1034 1.1384 
1.1805 1.22€0 1.2817 1.3424 1.4105 1.4871 
1.5731 1.6684 1.7736 1.8918 2.20229 2-1704 
2.3321 225169 2-/7161 2-2-9479 3-2-2018 324899 
3-8C88 4.1715 4.5826 5.0408 5-2-5616 6.1512 
6.8201 725845 6.4517 924439 10.5785 11.8784 


13.3689 15.0856 17.0596 19.3530 22.0055 25.0761 
28.6529 32.8206 37.6868 43.3807 50.0432 57.8872 
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Table 3. 
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Position and speed of ship for each record. 
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Heading, 
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Table 3. (cont. ) 
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Table 3. 
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(cont. ) 
Heading 
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Table 3. (cont. ) 


Position 
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Table 4. Data on short records for which less than 
600 points were available. 
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